Fast coarsening in unstable epitaxy with desorption 
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Homoepitaxial growth is unstable towards the formation of pyramidal mounds when interlayer 
transport is reduced due to activation barriers to hopping at step edges. Simulations of a lattice 
model and a continuum equation show that a small amount of desorption dramatically speeds up 
the coarsening of the mound array, leading to coarsening exponents between 1/3 and 1/2. The 
underlying mechanism is the faster growth of larger mounds due to their lower evaporation rate. 
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Growth of perfectly smooth epitaxial layers, e.g., 
for the fabrication of two-dimensional electron-gas het- 
erostructures, requires suppressing growth instabilities 
that lead to surface roughening. However, one can also 
try to harness the instabilities and make them produce 
regular arrays of tiny objects on the crystalline substrate, 
such as quantum wires and quantum dots. Much atten- 
tion has been paid to nanostructures created during het- 
eroepitaxial growth. It seems tempting to use also other 
types of instabilities associated with homoepitaxy. 

The most prominent of these kinetic instabilities is the 
creation of pyramidal features on the surface as a re- 
sult of additional activation barriers to hopping at step 
edges reducing interlayer transport [Qj2|. Theoretical 
|^-^| and experimental ^,|| research has shown that the 
array of pyramidal mounds coarsens in time with the av- 
erage mound size increasing according to a power 
law, £ cx i 1 / 2 . Under the standard assumption of negligi- 
ble desorption, which implies volume conservation for the 
growing film jl), the exponent 1/z was demonstrated to 
have an upper bound of 1/4 [ ^fj|Jllj| . Here we show, using 
kinetic Monte Carlo (KMC) simulations and numerical 
integration of a continuum equation of motion, that even 
a minute amount of desorption drastically changes this 
situation with 1/z increasing towards 1/2. Detailed in- 
vestigation of growth kinetics reveals that the reason is 
a dependence of the evaporation rate on the mound size, 
leading to faster growth of big mounds at the expense of 
small ones. 

Monte- Carlo simulations. To study epitaxial growth 
with desorption, we used a solid-on-solid KMC model in 
which the crystal is assumed to have a simple cubic struc- 
ture with neither bulk vacancies nor overhangs allowed. 
The basic processes included in our model are the depo- 
sition of atoms onto the surface at a rate F, their surface 
diffusion, and evaporation from the surface. The diffu- 
sion of surface adatoms is modeled as a nearest-neighbor 
hopping process at the rate kri = /coexp(— En/ksT), 
where Ed is the hopping barrier, T is the substrate 
temperature, and ks is Boltzmann's constant. The pre- 
factor ko is the attempt frequency of a surface adatom 



and is assigned the value 10 s . The barrier to hop- 
ping is given by E D = E s +nEN + (rii—nf)Q(ni-nf)E B 
where Es, En, and Eb are model parameters, n is the 
number of in-plane nearest neighbors before the hop, 
rii and n/ are the number of next-nearest neighbors in 
the planes beneath and above the hopping atom be- 
fore («i) and after (n/) a hop, and @{x)=l if x > 0, 
and otherwise (cf. Ref. Q for a more detailed de- 
scription of the model). The evaporation of a surface 
adatom occurs at the rate k cv — fco exp(— E cv /ksT), 
where E cv — Eq + nE^ with Eq being the energy for 
evaporation of a free surface adatom. 
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FIG. 1. Lateral mound size £(t) evolution for different 
model parameters and growth conditions obtained in KMC 
simulations. The model parameters used were i?o=1.9 eV 
at T=750 K, fij=0 (stars) and ,R,=3 (circles), £ =2.0 eV 
at T=750 K, Rt=3 (crosses), and E =1.75 eV at T=670 
K, Ri~3 (triangles). The slopes indicated are results of 
least-squares fitting. The error bars of the exponents (esti- 
mated from run-to-run variations) are of the order of 0.01. 
Inset: Size-dependent contribution to the evaporation rate of 
a single mound determined on lattices of size 21, 23, 27, 31, 
35, 41, 51, 61, 71, 81, and 101 with periodic boundary con- 
ditions. Data show average over 25 runs, in each of which 
about 1000 ML were deposited. The full line has the form 
Av = A/^ 1 ' 5 , and is consistent with coarsening exponents for 
the same model parameters. 
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The simulations were carried out on square 300x300 
to 600x600 lattices with periodic boundary conditions. 
The basic set of model parameters and growth conditions 
used was: E s = 1.54 eV, E N = 0.23 eV, E B = 0.175 eV, 
and F=l/6 monolayer (ML)/s (set I of Ref. |J). Un- 
der these conditions the equilibrium evaporation flux 
F oq = fc cxp[-(^ + 2E N )/k B T)] » (10- 3 - 10- 2 ) x F, 
and the actual desorption rate is a few times larger. The 
robustness of the observed behavior was tested by using 
different temperatures and evaporation barriers Eq, and 
by including the "incorporation radius" effect whereby 
the incoming atom is placed at the site with the highest 
number of lateral nearest neighbors within a square area 
of size Ri centered on the site of incidence [ 4J. 

Simulation results are shown in Fig. @. As the 
desorption rate increases, coarsening becomes faster and 
the coarsening exponent l/z becomes much bigger than 
0.19-0.26, the range of values observed in previous simu- 
lation work using the same model without desorption [Q . 
Even a rather small amount of desorption fl3|| thus dras- 
tically affects the coarsening exponent regardless of the 
details of the simulation model (such as the model pa- 
rameters and the incorporation radius). In the regions of 
fits shown in Fig. [j] the mound slope stays approximately 
constant, indicating that the asymptotic regime has been 
reached. An interesting feature of Fig. [j] is the crossover 
observed for the case Ri — after approx. 1000 ML 
were deposited. We discuss the underlying change in the 
mechanism of coarsening below. 




FIG. 2. Surface morphology in KMC simulations after ap- 
prox. 2000 ML have been deposited. The displayed part of the 
lattice is 300 x 300 (plane view, left) and f 00 x 100 (perspective 
plot, right). 

Fig. |^ shows plane and perspective views of the sur- 
face morphology after approx. 2000 ML have been de- 
posited. Pyramidal mounds are separated by narrow, 
deep troughs. The surface profile is clearly asymmetric 
with flat, rounded mound tops and sharp, deep valleys. 
For the case of Ri = 3, mounds are much shallower and 
bigger, having also more regular structure. In both cases, 
however, very fast coarsening is observed. 
Continuum equations. Further confirmation of the dra- 
matic effect of evaporation comes from continuum the- 
ory, where the surface is modeled by a smooth, space 
and time-dependent height function i?(x, £). Our start- 



ing point is the standard continuum equation for unstable 
epitaxy |7],[l(|, to which the leading order effect of des- 
orption is added in terms of a slope-dependent growth 
rate V(|VF|): 

°— = ~KA 2 H - V-[/(|Vtf | 2 )Vtf] + V(\VH\). (f) 

In the absence of desorption V(u) = F, the external flux. 
For a vicinal surface with step spacing £ (tilt u = a /£, 
with lattice constant a), the growth rate according to 
BCF theory reads Q VbcfM = (Fx s /£) t&nh(£/x s ) 
where x s = 1\J Dt is the desorption length, depending 
on the diffusion coefficient D and the desorption rate 1/r 
from a flat surface. 

To use the BCF-expression also for near singular sur- 
faces we introduce an effective, tilt-dependent step 
spacing £ c f[, which equals £ in the step flow regime, 
u 3> a j 'Id, and reduces to the terrace size or island dis- 
tance jf6| £n for small u. Desorption is considered a 
small, perturbative effect in the sense that 



£ D /x s < f , 



(2) 



which means that it is much more likely for an atom 
to be captured at a step than to desorb. Under this 
condition the terrace size £u should not be influenced 
by desorption. A plausible formula for £ e g is £ c g(u) 
= £d[1 + u 2 (£ D /a) 2 ]~ 1 /' 2 , and the growth rate is then 
V{u) = VbcF(4ff(«)). Because of (|), £ eS (u) < £ D < x s 
for all slopes, and therefore we can expand Vbcf to ob- 
tain 



V{u) « F[l - (a 2 /3)(l + u 2 (£ D /a) 2 Y 



(3) 

(l/3)« 2 ] 



The growth rate varies between V(0) = F[l 
for the singular surface and F for u 3> a/£o- 

In the second term on the right hand side of Eq. (|f|) 
we set ,f(u 2 ) = /q[1 — (u/mo) 2 ], which leads to a stable 
selected slope too a s is observed in our lattice sim- 
ulations ]17| |. We subtract the deposited film thickness, 
H — ► H — Ft and rescale time, lateral space, and height 
variables to arrive at the dimensionless form [fo| 



dh 



-A 2 h-V- 



l-(Vh) 2 )Vh 



a 2 /3 



1 + (V/>) 2 



(4) 



The one-dimensional version of (jl|) with an evaporation 
rate ~ (Vft.) 2 was considered in a different context by 
Emmott and Bray fll8[ . 

We integrated Eq. (|f|) numerically (for the method and 
system sizes see Jl9{) and found similar behavior as in 
the lattice model. After an initial fast increase of the 
lateral mound size £, the pattern coarsens as in the case 
without evaporation, £ ~ w ~ i 1 / 4 Jl0| (w is the mean 
square width of the surface, w 2 — (hr) where h — h—(h) 
denotes the height profile relative to its mean). This be- 
havior is transient and eventually crosses over to a fast 
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asymptotic increase of the mound size and the surface 
width as £ ~ w ~ t 1 / 2 . The mound size £(t) is shown 
in Figure || for values of a 2 /3 ranging from 10" 1 ' 2 to 
10 -3 / 2 , decreased by a factor 1CP 1 / 8 between succeed- 
ing curves. The transient t 1 ' 4 — regime is absent for the 
strongest evaporation (a 2 /3 = 10" 1 / 2 ) and becomes more 
pronounced as a is decreased. A similar crossover is ob- 
served in the KMC simulations with Ri = (stars in Fig. 
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FIG. 3. Lateral mound size £(£) for nine different evap- 
oration strengths « 2 /3 = l(r 1/2 , 1(T 5/8 , 10" 3/4 , 1(T 3/2 
obtained by numerical integration of the continuum equation. 
The transient regime £~i 1//4 persists until evaporation domi- 
nates the surface driven coarsening at t a ~a~ 4 , the crossover 
time to asymptotic fast coarsening £ ~ t 1 ^ 2 . Scaling plot in 
the inset shows time ixa 4 and length £ x a. 

The evaporation term in Eq. (Q) breaks the up-down 
symmetry (/i <-»• — h) pofl . When it is dominant (in the 
fast coarsening regime at late times) the surface mor- 
phology is asymmetric. The profile shown in Figure |] 
consists of conical mounds separated by well defined, nar- 
row valleys. Notice that there are no "negative mounds" . 
The greyscale plot shows the cellular arrangement of the 
cones. The observed features are very similar to results 
of simulations in Fig. ||. 



The origin of fast coarsening. Allowing for evaporation 
fundamentally changes the nature of the growth insta- 
bility, because it introduces a coupling between the lo- 
cal growth rate and the surface morphology. In particu- 
lar, one expects generically a dependence of the average 
growth rate of a mound on its size, while for conserved 
growth only the fluctuations in the growth rate are size 




FIG. 4. Surface profile from continuum equation (evapo- 
ration « 2 /3=10 -1/ ' 2 ) at late times when £ ~ t 1 ^ 2 . Coni- 
cal mounds (right) form a cellular structure visible in the 
greyscale representation (left). 



dependent |21|. To see how this can affect the coarsening 
law, assume that the growth velocity of a mound depends 
on its size £ as 



u(£) = vq — Av(£) where Au(£) ~ £ 



(5) 



with a positive prefactor, so that large mounds grow 
faster. If £ is the only macroscopic length in the system, 
the size differences between mounds are also of the order 
of £. The time scale on which a small mound is eliminated 
by its larger neighbors is then given by t% ~ w/Av(^), 
since the surface width w equals the typical height of 
mounds. Using that the mounds have a constant slope, 
to ~ £, it follows that £ ~ or z = 1 + v. Pro- 

vided v < 3 this violates the bound z > 4 obtained in 



the conserved case |10 



A well-known mechanism for a size-dependent growth 
rate is the Gibbs-Thomson-effect: For spherical droplets 
in equilibrium the evaporation rate is proportional to 
the curvature ~ l/£, hence within the Wilson- Frenkel- 
approximation (2^] the growth rate is of the form (||) with 
v = 1. The mounds in our lattice simulations are more 
conical in shape, with rather straight sides and rounded 
regions of lateral extent s» £o at the tips. Assuming that 
desorption occurs preferentially from the tip regions, the 
evaporation rate of a mound of size £ has a contribu- 
tion proportional to the ratio of the tip area ~ £ 2 D to 
the mound area ~ £ 2 , leading to v — 2 and 1/z = 1/3. 
For a more quantitative estimate, Au(£) was determined 
in a sequence of simulations on small square lattices of 
lateral size 21, 23, . . . up to 101 Q. As initial config- 
uration on each of them a single mound was prepared. 
It persisted during deposition of 1000 monolayers, and 
the average evaporation rate was determined from a se- 
quence of 25 runs. The data presented in Figure [l] show 
that Au(£) > and v = 1.5 ± 0.1, which is consistent 
with direct observations of coarsening on large lattices 
(cf. Fig. %. 

An analytical evaluation of At?(£) is possible for the 
continuum equation, which will also allow us to derive 
the scaling of the crossover times (cf. Fig. |3|) with a. 
We recall the surface profile of Figure |^. The cones show 
two lateral lengthscales: (i) their size £, which for late 
times is much larger than (ii) £d, the diameter of the 
tips and the valleys (=0(1) in our rescaled units) which 
is independent of £. 

Thus for a mound on a d-dimensional surface the frac- 
tion of the surface covered by the tip is (l/£) d . The 
surrounding trough has codimension one and a relative 
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weight l/£, while the major part of the surface consists 
of the sloped sides of the conical mounds. 

Evaporation is less pronounced on the mounds' sides 
whereas it is enhanced by an amount of order a 2 on the 
horizontal parts, i.e. on the tip and in the surrounding 
valley. As a consequence equation (||) also holds for the 
continuum equation, where vq is the evaporation rate 
on the mounds' sides, and the enhanced mass loss from 
small mounds is mainly due to the surrounding valley, 
i.e. Au(£) ~ a 2 /£- So the timescale for mound coales- 
cence is ~ w/Av(£) ~ £, 2 /a 2 (due to the stable slope, 
and it follows that £ ~ at 1 / 2 . Incidentally, the 
same coarsening law was found in the one-dimensional 
case @. 

The initial increase £ ~ t 1 / 4 is not due to evaporation 
and thus the same for all values of a (see Figure |J). To- 
gether with the late time behavior £ ~ at 1 / 2 this yields 
the estimate t a ^a~ 4 — (ie s /£d) 4 for the time at which 
evaporation begins to dominate the coarsening process. 
Rescaling time as t/t a and length as £/t^/ 4 , and omitting 
the initial fast increase puts all curves and in particular 
the crossover times t a on top of each other, as shown in 
the inset of Figure 0. 

We emphasize that fast coarsening for the continuum 
equation is due to the dominance of evaporation from 
valleys compared to the rest of the mounds. Direct in- 
spection shows that in the KMC simulations more atoms 
in fact evaporate from the upper parts of the mounds. 
This explains why the coarsening exponent observed for 
the lattice model is smaller than 1/2: Given enhanced 
evaporation only on the tips, Au(£) is of the order of 
a 2 /£ d , leading to £ ~ (a 2 t) 1 ^ d+1 ^ in d dimensions, hence 
z = 3 for d = 2 as argued previously. To improve on this 
estimate, more detailed information about the shape of 
mounds and its coupling to the evaporation rate would 
be needed. It is nevertheless interesting to note that 
the coarsening exponent 1/z = 1) is always larger 

than the value 1/z = l/(d+2) obtained for noise-induced 
coarsening pl[|, indicating that our conclusions will not 
be modified by shot noise. 

In summary, we have identified a general mechanism 
for fast mound coarsening in unstable growth with des- 
orption. While the detailed appearance of the effect is 
different in the lattice model as compared to the con- 
tinuum equation, in both cases the key feature is the 
dependence of growth rate on mound size. This gives 
us confidence that the phenomenon is robust and will be 
observed under suitable experimental conditions. 
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